Genetic determinants of blood pressure and heart rate identified through ENU-induced mutagenesis with automated meiotic mapping.

We used N-ethyl-N-nitrosurea-induced germline mutagenesis combined with automated meiotic mapping to identify specific systolic blood pressure (SBP) and heart rate (HR) determinant loci. We analyzed 43,627 third-generation (G3) mice from 841 pedigrees to assess the effects of 45,378 variant alleles within 15,760 genes, in both heterozygous and homozygous states. We comprehensively tested 23% of all protein-encoding autosomal genes and found 87 SBP and 144 HR (with 7 affecting both) candidates exhibiting detectable hypomorphic characteristics. Unexpectedly, only 18 of the 87 SBP genes were previously known, while 26 of the 144 genes linked to HR were previously identified. Furthermore, we confirmed the influence of two genes on SBP regulation and three genes on HR control through reverse genetics. This underscores the importance of our research in uncovering genes associated with these critical cardiovascular risk factors and illustrate the effectiveness of germline mutagenesis for defining key determinants of polygenic phenotypes that must be studied in an intact organism.

Obesity caused by an OVOL2 mutation reveals dual roles of OVOL2 in promoting thermogenesis and limiting white adipogenesis.

Using random germline mutagenesis in mice, we identified a viable hypomorphic allele (boh) of the transcription-factor-encoding gene Ovol2 that resulted in obesity, which initially developed with normal food intake and physical activity but decreased energy expenditure. Fat weight was dramatically increased, while lean weight was reduced in 12-week-old boh homozygous mice, culminating by 24 weeks in massive obesity, hepatosteatosis, insulin resistance, and diabetes. The Ovol2boh/boh genotype augmented obesity in Lepob/ob mice, and pair-feeding failed to normalize obesity in Ovol2boh/boh mice. OVOL2-deficient mice were extremely cold intolerant. OVOL2 is essential for brown/beige adipose tissue-mediated thermogenesis. In white adipose tissues, OVOL2 limited adipogenesis by blocking C/EBPα engagement of its transcriptional targets. Overexpression of OVOL2 in adipocytes of mice fed with a high-fat diet reduced total body and liver fat and improved insulin sensitivity. Our data reveal that OVOL2 plays dual functions in thermogenesis and adipogenesis to maintain energy balance.

Modulation of autoimmune diabetes by N-ethyl-N-nitrosourea- induced mutations in non-obese diabetic mice.

Genetic association studies of type 1 diabetes (T1D) in humans, and in congenic non-obese diabetic (NOD) mice harboring DNA segments from T1D-resistant mice, face the challenge of assigning causation to specific gene variants among many within loci that affect disease risk. Here, we created random germline mutations in NOD/NckH mice and used automated meiotic mapping to identify mutations modifying T1D incidence and age of onset. In contrast with association studies in humans or congenic NOD mice, we analyzed a relatively small number of genetic changes in each pedigree, permitting implication of specific mutations as causative. Among 844 mice from 14 pedigrees bearing 594 coding/splicing changes, we identified seven mutations that accelerated T1D development, and five that delayed or suppressed T1D. Eleven mutations affected genes not previously known to influence T1D (Xpnpep1, Herc1, Srrm2, Rapgef1, Ppl, Zfp583, Aldh1l1, Col6a1, Ccdc13, Cd200r1, Atrnl1). A suppressor mutation in Coro1a validated the screen. Mutagenesis coupled with automated meiotic mapping can detect genes in which allelic variation influences T1D susceptibility in NOD mice. Variation of some of the orthologous/paralogous genes may influence T1D susceptibility in humans.

RNPS1 inhibits excessive tumor necrosis factor/tumor necrosis factor receptor signaling to support hematopoiesis in mice.

Null mutations of spliceosome components or cofactors are homozygous lethal in eukaryotes, but viable hypomorphic mutations provide an opportunity to understand the physiological impact of individual splicing proteins. We describe a viable missense allele (F181I) of Rnps1 encoding an essential regulator of splicing and nonsense-mediated decay (NMD), identified in a mouse genetic screen for altered immune cell development. Homozygous mice displayed a stem cell­intrinsic defect in hematopoiesis of all lineages due to excessive apoptosis induced by tumor necrosis factor (TNF)­dependent death signaling. Numerous transcript splice variants containing retained introns and skipped exons were detected at elevated frequencies in Rnps1F181I/F181I splenic CD8+ T cells and hematopoietic stem cells (HSCs), but NMD appeared normal. Strikingly, Tnf knockout rescued all hematopoietic cells to normal or near-normal levels in Rnps1F181I/F181I mice and dramatically reduced intron retention in Rnps1F181I/F181I CD8+ T cells and HSCs. Thus, RNPS1 is necessary for accurate splicing, without which disinhibited TNF signaling triggers hematopoietic cell death.

De novo germline mutation in the dual specificity phosphatase 10 gene accelerates autoimmune diabetes.

Insulin-dependent or type 1 diabetes (T1D) is a polygenic autoimmune disease. In humans, more than 60 loci carrying common variants that confer disease susceptibility have been identified by genome-wide association studies, with a low individual risk contribution for most variants excepting those of the major histocompatibility complex (MHC) region (40 to 50% of risk); hence the importance of missing heritability due in part to rare variants. Nonobese diabetic (NOD) mice recapitulate major features of the human disease including genetic aspects with a key role for the MHC haplotype and a series of Idd loci. Here we mapped in NOD mice rare variants arising from genetic drift and significantly impacting disease risk. To that aim we established by selective breeding two sublines of NOD mice from our inbred NOD/Nck colony exhibiting a significant difference in T1D incidence. Whole-genome sequencing of high (H)- and low (L)-incidence sublines (NOD/NckH and NOD/NckL) revealed a limited number of subline-specific variants. Treating age of diabetes onset as a quantitative trait in automated meiotic mapping (AMM), enhanced susceptibility in NOD/NckH mice was unambiguously attributed to a recessive missense mutation of Dusp10, which encodes a dual specificity phosphatase. The causative effect of the mutation was verified by targeting Dusp10 with CRISPR-Cas9 in NOD/NckL mice, a manipulation that significantly increased disease incidence. The Dusp10 mutation resulted in islet cell down-regulation of type I interferon signature genes, which may exert protective effects against autoimmune aggression. De novo mutations akin to rare human susceptibility variants can alter the T1D phenotype.

Thousands of induced germline mutations affecting immune cells identified by automated meiotic mapping coupled with machine learning.

Forward genetic studies use meiotic mapping to adduce evidence that a particular mutation, normally induced by a germline mutagen, is causative of a particular phenotype. Particularly in small pedigrees, cosegregation of multiple mutations, occasional unawareness of mutations, and paucity of homozygotes may lead to erroneous declarations of cause and effect. We sought to improve the identification of mutations causing immune phenotypes in mice by creating Candidate Explorer (CE), a machine-learning software program that integrates 67 features of genetic mapping data into a single numeric score, mathematically convertible to the probability of verification of any putative mutation-phenotype association. At this time, CE has evaluated putative mutation-phenotype associations arising from screening damaging mutations in ∼55% of mouse genes for effects on flow cytometry measurements of immune cells in the blood. CE has therefore identified more than half of genes within which mutations can be causative of flow cytometric phenovariation in Mus musculus The majority of these genes were not previously known to support immune function or homeostasis. Mouse geneticists will find CE data informative in identifying causative mutations within quantitative trait loci, while clinical geneticists may use CE to help connect causative variants with rare heritable diseases of immunity, even in the absence of linkage information. CE displays integrated mutation, phenotype, and linkage data, and is freely available for query online.

SLFN2 protection of tRNAs from stress-induced cleavage is essential for T cell-mediated immunity.

Reactive oxygen species (ROS) increase in activated T cells because of metabolic activity induced to support T cell proliferation and differentiation. We show that these ROS trigger an oxidative stress response that leads to translation repression. This response is countered by Schlafen 2 (SLFN2), which directly binds transfer RNAs (tRNAs) to protect them from cleavage by the ribonuclease angiogenin. T cell-specific SLFN2 deficiency results in the accumulation of tRNA fragments, which inhibit translation and promote stress-granule formation. Interleukin-2 receptor ß (IL-2Rß) and IL-2Rγ fail to be translationally up-regulated after T cell receptor stimulation, rendering SLFN2-deficient T cells insensitive to interleukin-2's mitogenic effects. SLFN2 confers resistance against the ROS-mediated translation-inhibitory effects of oxidative stress normally induced by T cell activation, permitting the robust protein synthesis necessary for T cell expansion and immunity.

N4BP1 negatively regulates NF-κB by binding and inhibiting NEMO oligomerization.

Many immune responses depend upon activation of NF-κB, an important transcription factor in the elicitation of a cytokine response. Here we show that N4BP1 inhibits TLR-dependent activation of NF-κB by interacting with the NF-κB signaling essential modulator (NEMO, also known as IκB kinase γ) to attenuate NEMO-NEMO dimerization or oligomerization. The UBA-like (ubiquitin associated-like) and CUE-like (ubiquitin conjugation to ER degradation-like) domains in N4BP1 mediate interaction with the NEMO COZI domain. Both in vitro and in mice, N4bp1 deficiency specifically enhances TRIF-independent (TLR2, TLR7, or TLR9-mediated) but not TRIF-dependent (TLR3 or TLR4-mediated) NF-κB activation, leading to increased production of proinflammatory cytokines. In response to TLR4 or TLR3 activation, TRIF causes activation of caspase-8, which cleaves N4BP1 distal to residues D424 and D490 and abolishes its inhibitory effect. N4bp1-/- mice also have diminished numbers of T cells in the peripheral blood. Our work identifies N4BP1 as an inhibitory checkpoint protein that must be overcome to activate NF-κB, and a TRIF-initiated caspase-8-dependent mechanism by which this is accomplished.

Dominant atopy risk mutations identified by mouse forward genetic analysis.

BACKGROUND: Atopy, the overall tendency to become sensitized to an allergen, is heritable but seldom ascribed to mutations within specific genes. Atopic individuals develop abnormally elevated IgE responses to immunization with potential allergens. To gain insight into the genetic causes of atopy, we carried out a forward genetic screen for atopy in mice. METHODS: We screened mice carrying homozygous and heterozygous N-ethyl-N-nitrosourea (ENU)-induced germline mutations for aberrant antigen-specific IgE and IgG1 production in response to immunization with the model allergen papain. Candidate genes were validated by independent gene mutation. RESULTS: Of 31 candidate genes selected for investigation, the effects of mutations in 23 genes on papain-specific IgE or IgG1 were verified. Among the 20 verified genes influencing the IgE response, eight were necessary for the response, while 12 repressed IgE. Nine genes were not previously implicated in the IgE response. Fifteen genes encoded proteins contributing to IgE class switch recombination or B-cell receptor signaling. The precise roles of the five remaining genes (Flcn, Map1lc3b, Me2, Prkd2, and Scarb2) remain to be determined. Loss-of-function mutations in nine of the 12 genes limiting the IgE response were dominant or semi-dominant for the IgE phenotype but did not cause immunodeficiency in the heterozygous state. Using damaging allele frequencies for the corresponding human genes and in silico simulations (Monte Carlo) of undiscovered atopy mutations, we estimated the percentage of humans with heterozygous atopy risk mutations. CONCLUSIONS: Up to 37% of individuals may be heterozygous carriers for at least one dominant atopy risk mutation.

Mutual inhibition between Prkd2 and Bcl6 controls T follicular helper cell differentiation.

T follicular helper cells (TFH) participate in germinal center (GC) development and are necessary for B cell production of high-affinity, isotype-switched antibodies. In a forward genetic screen, we identified a missense mutation in Prkd2, encoding the serine/threonine kinase protein kinase D2, which caused elevated titers of immunoglobulin E (IgE) in the serum. Subsequent analysis of serum antibodies in mice with a targeted null mutation of Prkd2 demonstrated polyclonal hypergammaglobulinemia of IgE, IgG1, and IgA isotypes, which was exacerbated by the T cell-dependent humoral response to immunization. GC formation and GC B cells were increased in Prkd2-/- spleens. These effects were the result of excessive cell-autonomous TFH development caused by unrestricted Bcl6 nuclear translocation in Prkd2-/- CD4+ T cells. Prkd2 directly binds to Bcl6, and Prkd2-dependent phosphorylation of Bcl6 is necessary to constrain Bcl6 to the cytoplasm, thereby limiting TFH development. In response to immunization, Bcl6 repressed Prkd2 expression in CD4+ T cells, thereby committing them to TFH development. Thus, Prkd2 and Bcl6 form a mutually inhibitory positive feedback loop that controls the stable transition from naïve CD4+ T cells to TFH during the adaptive immune response.

Syndromic immune disorder caused by a viable hypomorphic allele of spliceosome component Snrnp40.

We report a new immunodeficiency disorder in mice caused by a viable hypomorphic mutation of Snrnp40, an essential gene encoding a subunit of the U5 small nuclear ribonucleoprotein (snRNP) complex of the spliceosome. Snrnp40 is ubiquitous but strongly expressed in lymphoid tissue. Homozygous mutant mice showed hypersusceptibility to infection by murine cytomegalovirus and multiple defects of lymphoid development, stability and function. Cell-intrinsic defects of hematopoietic stem cell differentiation also affected homozygous mutants. SNRNP40 deficiency in primary hematopoietic stem cells or T cells or the EL4 cell line increased the frequency of splicing errors, mostly intron retention, in several hundred messenger RNAs. Altered expression of proteins associated with immune cell function was also observed in Snrnp40-mutant cells. The immunological consequences of SNRNP40 deficiency presumably result from cumulative, moderate effects on processing of many different mRNA molecules and secondary reductions in the expression of critical immune proteins, yielding a syndromic immune disorder.

Adjuvant effect of the novel TLR1/TLR2 agonist Diprovocim synergizes with anti-PD-L1 to eliminate melanoma in mice.

Successful cancer immunotherapy entails activation of innate immune receptors to promote dendritic cell (DC) maturation, antigen presentation, up-regulation of costimulatory molecules, and cytokine secretion, leading to activation of tumor antigen-specific cytotoxic T lymphocytes (CTLs). Here we screened a synthetic library of 100,000 compounds for innate immune activators using TNF production by THP-1 cells as a readout. We identified and optimized a potent human and mouse Toll-like receptor (TLR)1/TLR2 agonist, Diprovocim, which exhibited an EC50 of 110 pM in human THP-1 cells and 1.3 nM in primary mouse peritoneal macrophages. In mice, Diprovocim-adjuvanted ovalbumin immunization promoted antigen-specific humoral and CTL responses and synergized with anti-PD-L1 treatment to inhibit tumor growth, generating long-term antitumor memory, curing or prolonging survival of mice engrafted with the murine melanoma B16-OVA. Diprovocim induced greater frequencies of tumor-infiltrating leukocytes than alum, of which CD8 T cells were necessary for the antitumor effect of immunization plus anti-PD-L1 treatment.

Probability of phenotypically detectable protein damage by ENU-induced mutations in the Mutagenetix database.

Computational inference of mutation effects is necessary for genetic studies in which many mutations must be considered as etiologic candidates. Programs such as PolyPhen-2 predict the relative severity of damage caused by missense mutations, but not the actual probability that a mutation will reduce/eliminate protein function. Based on genotype and phenotype data for 116,330 ENU-induced mutations in the Mutagenetix database, we calculate that putative null mutations, and PolyPhen-2-classified "probably damaging", "possibly damaging", or "probably benign" mutations have, respectively, 61%, 17%, 9.8%, and 4.5% probabilities of causing phenotypically detectable damage in the homozygous state. We use these probabilities in the estimation of genome saturation and the probability that individual proteins have been adequately tested for function in specific genetic screens. We estimate the proportion of essential autosomal genes in Mus musculus (C57BL/6J) and show that viable mutations in essential genes are more likely to induce phenotype than mutations in non-essential genes.

HCFC2 is needed for IRF1- and IRF2-dependent Tlr3 transcription and for survival during viral infections.

Transcriptional regulation of numerous interferon-regulated genes, including Toll-like receptor 3 (Tlr3), which encodes an innate immune sensor of viral double-stranded RNA, depends on the interferon regulatory factor 1 (IRF1) and IRF2 transcription factors. We detected specific abrogation of macrophage responses to polyinosinic-polycytidylic acid (poly(I:C)) resulting from three independent N-ethyl-N-nitrosourea-induced mutations in host cell factor C2 (Hcfc2). Hcfc2 mutations compromised survival during influenza virus and herpes simplex virus 1 infections. HCFC2 promoted the binding of IRF1 and IRF2 to the Tlr3 promoter, without which inflammatory cytokine and type I IFN responses to the double-stranded RNA analogue poly(I:C) are reduced in mouse macrophages. HCFC2 was also necessary for the transcription of a large subset of other IRF2-dependent interferon-regulated genes. Deleterious mutations of Hcfc2 may therefore increase susceptibility to diverse infectious diseases.

Skin-specific regulation of SREBP processing and lipid biosynthesis by glycerol kinase 5.

The recessive N-ethyl-N-nitrosourea-induced phenotype toku is characterized by delayed hair growth, progressive hair loss, and excessive accumulation of dermal cholesterol, triglycerides, and ceramides. The toku phenotype was attributed to a null allele of Gk5, encoding glycerol kinase 5 (GK5), a skin-specific kinase expressed predominantly in sebaceous glands. GK5 formed a complex with the sterol regulatory element-binding proteins (SREBPs) through their C-terminal regulatory domains, inhibiting SREBP processing and activation. In Gk5toku/toku mice, transcriptionally active SREBPs accumulated in the skin, but not in the liver; they were localized to the nucleus and led to elevated lipid synthesis and subsequent hair growth defects. Similar defective hair growth was observed in kinase-inactive GK5 mutant mice. Hair growth defects of homozygous toku mice were partially rescued by treatment with the HMG-CoA reductase inhibitor simvastatin. GK5 exists as part of a skin-specific regulatory mechanism for cholesterol biosynthesis, independent of cholesterol regulation elsewhere in the body.

IgD class switching is initiated by microbiota and limited to mucosa-associated lymphoid tissue in mice.

Class-switch recombination (CSR) alters the Ig isotype to diversify antibody effector functions. IgD CSR is a rare event, and its regulation is poorly understood. We report that deficiency of 53BP1, a DNA damage-response protein, caused age-dependent overproduction of secreted IgD resulting from increased IgD CSR exclusively within B cells of mucosa-associated lymphoid tissues. IgD overproduction was dependent on activation-induced cytidine deaminase, hematopoietic MyD88 expression, and an intact microbiome, against which circulating IgD, but not IgM, was reactive. IgD CSR occurred via both alternative nonhomologous end-joining and homologous recombination pathways. Microbiota-dependent IgD CSR also was detected in nasal-associated lymphoid tissue of WT mice. These results identify a pathway, present in WT mice and hyperactivated in 53BP1-deficient mice, by which microbiota signal via Toll-like receptors to elicit IgD CSR.

NLRP3 activation and mitosis are mutually exclusive events coordinated by NEK7, a new inflammasome component.

The NLRP3 inflammasome responds to microbes and danger signals by processing and activating proinflammatory cytokines, including interleukin 1ß (IL-1ß) and IL-18. We found here that activation of the NLRP3 inflammasome was restricted to interphase of the cell cycle by NEK7, a serine-threonine kinase previously linked to mitosis. Activation of the NLRP3 inflammasome required NEK7, which bound to the leucine-rich repeat domain of NLRP3 in a kinase-independent manner downstream of the induction of mitochondrial reactive oxygen species (ROS). This interaction was necessary for the formation of a complex containing NLRP3 and the adaptor ASC, oligomerization of ASC and activation of caspase-1. NEK7 promoted the NLRP3-dependent cellular inflammatory response to intraperitoneal challenge with monosodium urate and the development of experimental autoimmune encephalitis in mice. Our findings suggest that NEK7 serves as a cellular switch that enforces mutual exclusivity of the inflammasome response and cell division.

Real-time resolution of point mutations that cause phenovariance in mice.

With the wide availability of massively parallel sequencing technologies, genetic mapping has become the rate limiting step in mammalian forward genetics. Here we introduce a method for real-time identification of N-ethyl-N-nitrosourea-induced mutations that cause phenotypes in mice. All mutations are identified by whole exome G1 progenitor sequencing and their zygosity is established in G2/G3 mice before phenotypic assessment. Quantitative and qualitative traits, including lethal effects, in single or multiple combined pedigrees are then analyzed with Linkage Analyzer, a software program that detects significant linkage between individual mutations and aberrant phenotypic scores and presents processed data as Manhattan plots. As multiple alleles of genes are acquired through mutagenesis, pooled "superpedigrees" are created to analyze the effects. Our method is distinguished from conventional forward genetic methods because it permits (1) unbiased declaration of mappable phenotypes, including those that are incompletely penetrant (2), automated identification of causative mutations concurrent with phenotypic screening, without the need to outcross mutant mice to another strain and backcross them, and (3) exclusion of genes not involved in phenotypes of interest. We validated our approach and Linkage Analyzer for the identification of 47 mutations in 45 previously known genes causative for adaptive immune phenotypes; our analysis also implicated 474 genes not previously associated with immune function. The method described here permits forward genetic analysis in mice, limited only by the rates of mutant production and screening.